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Scaling arguments are developed for the load balance in hydrodynamic lubrication, and applied to 
non-Newtonian lubricants with a shear-thinning rheology typical of a structured liquid. It is argued 
that the shear thinning regime may be mechanically unstable in lubrication flow, and consequently 
the Stribeck (friction) curve should be discontinuous, with possible hysteresis. Further analysis 
suggests that normal stress and flow transience (stress overshoot) do not destroy this basic picture, 
although they may provide stabilising mechanisms at higher shear rates. Extensional viscosity is 
also expected to be insignificant unless the Trouton ratio is large. A possible application to recent 
theories of shear thickening in non-Brownian particulate suspensions is indicated. 
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I. INTRODUCTION 

Lubrication is usually considered to be a complex, 
multi-scale, mechanical engineering problem Q}[3], de¬ 
manding sophisticated numerical approaches for its so¬ 
lution @ 10 . Indeed, with few exceptions [6], pure hydro- 
dynamic lubrication is rarely studied; rather the focus is 
more usually on elasto-hydrodynamic lubrication (EHL) 
or boundary layer lubrication. Given this, the utility of 
a simplified scaling theory for pure hydrodynamic lubri¬ 
cation may be questioned. However there are numer¬ 
ous areas to which lubrication is additionally relevant, 
for example in understanding the origins of ‘mouth-feel’ 
and the sensory properties of foods [7], and in ‘psycho¬ 
rheology’ and the perception of skin care products [5]. 
These diverse applications differ from the traditional en¬ 
gineering problem in that the lubricants are often struc¬ 
tured liquids or soft solids with a non-Newtonian rheol¬ 
ogy, and the gaps are large enough that pure hydrody¬ 
namic lubrication may be relevant. This motivates the 
development of the present scaling theory, which aims to 
provide insights and guidance for these unconventional 
application areas, and perhaps even reveal qualitatively 
new phenomena. 

To begin with, let me outline the fundamental mecha¬ 
nism of pure hydrodynamic lubrication [5] [TO] . In a fully 
lubricated conjunction (Fig. [I]), mass conservation within 
the converging and diverging wedges induces a Poiseuille- 
like contribution to the entrainment flow, superimposed 
on a Couette-type shearing motion. There is a corre¬ 
sponding emergent pressure distribution (Fig. [I] inset), 
which supports the load. The ratio of the load to the 
lateral sliding force defines the friction coefficient fj.. 

As we shortly shall see, Reynolds lubrication theory 
predicts /r is proportional to the Sommerfeld number pQ 
S = r]UR/W, where rj is the lubricant viscosity, U is 
the sliding velocity, R is a length scale characterising the 
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curvature of the surfaces, and W is the normal load. The 
overall p(S) behaviour is often summarised in the semi- 
empirical Stribeck curve (Fig. [2]). Note that pure hydro- 
dynamic scaling breaks down as S —> 0 since the surfaces 
come more and more into close contact where EHL be¬ 
comes increasingly important. This is reflected in the be¬ 
haviour observed empirically in the Stribeck curve. The 
present scaling theory should be applied, and interpreted, 
within this wider context. 

Before embarking on the detailed development, one 
further general remark should be made. In a geometry 
which posesses reflection symmetry, Reynolds lubrication 
theory predicts the lubrication pressure should be skew- 
symmetric about the minimum gap (Fig. [ljinset). There¬ 
fore, the integrated lubrication pressure should vanish. 
In actuality it has long been recognised that some addi¬ 
tional physics intervenes to knock down this result. For 
instance in the trailing edge where the lubrication pres¬ 
sure is sub-ambient (shaded area in Fig. [I] inset), the 
free surface may separate HH, or the fluid may cavitate 
fl2l . These considerations make the exact solution to 
the problem dependent on the nature of the additional 



FIG. 1. The geometry of a sphere and a flat. It is supposed 
that the flat surface is sliding at a velocity U underneath 
the sphere. Inset: the skew-symmetric lubrication pressure 
distribution. 
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FIG. 2. Stribeck curve (friction coefficient as a function of 
Sommerfeld number): EHL is elasto-hydrodynamic lubrica¬ 
tion. In the hydrodynamic regime, the friction coefficient is 
proportional to the Sommerfeld number (dashed line). 


physics. However a simple and commonly used prescrip¬ 
tion is to discard the contribution from the sub-ambient 
pressure region (the so-called half-Sommerfeld boundary 
condition El)- I shall tacitly use this assumption below. 

Keeping all this in mind, in the next section I shall 
develop the scaling theory for Newtonian hydrodynamic 
lubrication, and in the following section I shall explore 
the implications for non-Newtonian liquids. 


II. SCALING THEORY: NEWTONIAN CASE 


The gap between the lubricated surfaces, as shown in 
Fig. [I] for example, can always be represented in the con¬ 
junction region by a parabolic profile 

h = + 2R ' (•*■) 

In this ho is the minimum gap, r is the radial distance 
from the minimum gap, and I? is a measure of the radius 
of curvature of the surfaces (e. g. the sphere radius in 
Fig. 0- The radius within which h remains of order ho 
defines a length scale r 0 ~ ( Rh 0 j 1 / 2 . This length scale 
sets the size of the conjunction region and is the key to 
the development of the scaling theory. 

As already alluded to, and described in more detail 
in the Appendix, mass conservation induces a Poiseuille- 
type flow of magnitude U superimposed on the Couette- 
type shearing motion. The matching pressure distribu¬ 
tion is such that U ~ {ho/rfj’Vp where Vp ~ p/r o is pres¬ 
sure gradient. Hence the Reynolds lubrication pressure 
has a magnitude 


P 


rjUro 

hi 


( 2 ) 


The integrated pressure (e. g. the unshaded area in the 
Fig. [l] inset) then corresponds to a normal force 


W o 


rjUrf 


h o 


^ ~ p[/ii 3/,2 /lg 1/2 . 


( 3 ) 



FIG. 3. Load curve in the Newtonian case. The solid line is 
the normal force generated by the lubrication pressure. The 
filled circle is a stable fixed point, under a given load. 


In this I have used ro ~ (Rho) 1 ^ 2 for the second step. 
This integrated pressure balances the load when pr^ ~ 
W, or alternatively, using Eq. ([3|, when 


ho 

~R 



( 4 ) 


More careful treatments restore a prefactor to this, for 
example Kapitza m derived ho = (727r 2 /25) RS 2 « 
28.4 RS 2 , and Hamrock [2] has ho ~ 34.8 RS 2 . 

A point not often stressed, but which will be critical 
in the sequel, is that the lubricated conjunction is me¬ 
chanically stable at the load balance condition. To see 
this suppose that the minimum gap ho has not yet taken 
its steady state value. The load W presses down on the 
conjunction, and is opposed by the lubrication pressure 
which exerts a force of order prl given by Eq. © , which 
I note is a decreasing function of ho- Therefore if the gap 
is too large, the lubrication pressure does not support 
the load and the gap closes. Conversely, if the gap is too 
small, the lubrication pressure overcompensates for the 
load and the gap opens up. This stabilising mechanism is 
shown schematically in Fig. [3j where the solid line is the 
normal force arising from the lubrication pressure. In the 
language of dynamical systems theory, the filled circle in 
this diagram is a stable fixed point. 

I turn now to the Stribeck curve. The tangential wall 
stress in the conjunction zone is of order pU /ho (see Ap¬ 
pendix). Multiplying by the area of the conjunction (rg) 
gives the tangential force T ~ r]Url/ho■ Hence the fric¬ 
tion coefficient 


= T_ rjUrl h_UR 

M W~ Who ~ W ^ ’ 


( 5 ) 


(using the geometrical condition rg ~ Rho). This ex¬ 
plains why /i (X S in the hydrodynamic part of the 
Stribeck curve. As already noted, the actual friction co¬ 
efficient does not decrease indefinitely as S —► 0 since the 
gap shrinks, i. e. h 0 ~ S 2 from Eq. Q. Eventually the 
gap becomes very small and the friction coefficient devi¬ 
ates from the pure hydrodynamic law at the point where 
one enters the EHL regime m 
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FIG. 4. (a) Load curves in the non-Newtonian strong shear 
thinning case. The solid line in is the normal force generated 
by the lubrication pressure. A filled (open) circle is a stable 
(unstable) fixed point, under a given load. The fixed points 
disappear in a saddle-node bifurcation as the load is increased, 
(b) The corresponding Stribeck curve is predicted to jump 
discontinuously at the critical load, perhaps into a boundary 
lubrication regime. 


III. NON-NEWTONIAN CASE 

I now propose the obvious extensions of the above scal¬ 
ing arguments to incorporate non-linear rheological ef¬ 
fects of interest for non-Newtonian lubricants. For the 
time being I shall leave aside the question of normal 
stress, flow transience, and extensional flow, and focus 
first on the typical shear thinning behaviour found in 
structured liquids. 


A. Shear thinning 


I shall represent the non-linear viscosity in this case by 
the simple model, cf. [5], 


Vo 

1 +(yr) a 


( 6 ) 


FIG. 5. (a) Lubrication pressure in the case of shear thinning 
with a second Newtonian plateau, (b) Consequential hystere¬ 
sis in the Stribeck curve. 


a »o.8-i nun, although I should note that a > 1 
would correspond to shear banding. Eq. ([ 6 | omits the 
second Newtonian plateau which is expected to obtain 
at very high shear rates; this will be added informally to 
the discussion at the end. 

Inserting Eq. ^ into Eq. ^ gives a revised estimate 
for the normal force produced by the lubrication pressure, 


P r o 


Vo 

1 + (jr) a 


x UR 3/2 h~ 1/2 


(7) 


where the shear rate is estimated by 7 ~ U/ho- The per¬ 
haps naive and rather drastic assumptions being made 
here will be reviewed later. Since the shear rate is a 
decreasing function of the minimum gap, there are two 
regimes. For large gap, the Weissenberg number is small 
(7T < 1) and one recovers the Newtonian behaviour seen 
already, with 77 = 770 being the viscosity in the first New¬ 
tonian plateau. 

For small gap, one enters the high Weissenberg number 
regime (qr > 1) in which case Eq. (Jt|) becomes 

Prl ~ 77o(7r)"“ UR 3 / 2 h~ 1/2 ~ /7“- 1/2 • (8) 


where 770 is the low shear viscosity (first Newtonian 
plateau), 7 is the shear rate, r is a characteristic relax¬ 
ation time, and a is an exponent characterising the rate 
of decrease of the viscosity in the shear-thinning regime. 
The dimensionless quantity qr is the Weissenberg number 
and is a measure of the extent to which the non-linear 
regime has been penetrated. For concreteness, a typi¬ 
cal value of the exponent for polymer melt rheology is 


Only the Tio-dependence has been retained in the final 
step. If a < 1/2 (weak shear thinning), then Eq. ([ 8 ]) is a 
decreasing function of ho and one would again expect the 
gap to be mechanically stable, albeit with a dependence 
on a modified or instantaneous Sommerfeld number. On 
the other hand, if a > 1/2 (strong shear thinning), pr% 
in Eq. ( 8 | is an increasing function of ho and therefore 
the load balance in this regime will be unstable. 
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I shall take this latter case (strong shear thinning) to be 
the one of most interest (for polymer melts, for instance). 
The two regimes in the load balance are illustrated in 
Fig. [§r. The left hand branch is the high Weissenberg 
number regime where shear thinning takes place. In this 
branch, the open circle corresponds to an unstable fixed 
point at load W ±. To the left of the open circle, the lubri¬ 
cation pressure is insufficient to support the load and the 
gap closes until interrupted by some additional physics. 
To the right of the open circle, the lubrication pressure 
forces the surfaces to move apart until one reaches the 
stable fixed point (filled circle) in the low Weissenberg 
number regime. Now consider increasing the load to Wi- 
In this situation, the lubrication pressure is never suffi¬ 
cient to keep the surfaces apart and the gap closes un¬ 
til one enters the mixed or boundary layer lubrication 
regime. Somewhere in between W\ and W 2 is a critical 
load at which the stable and unstable fixed points merge 
in a saddle-node bifurcation. 

The implication for the Stribeck curve is as follows. In 
the low Weissenberg number regime, the analysis goes 
through as for the Newtonian case, with 770 featuring as 
the viscosity. The Weissenberg number yr ~ S~ 2 , and 
increases as the Sommerfeld number shrinks. Mechanical 
stability is lost at the point where one enters the shear 
thinning regime (i. e. at the above saddle-node bifurca¬ 
tion). At this point the surfaces should jump into (near) 
contact. Correspondingly the Stribeck curve should show 
a discontinuous jump, most likely into a boundary lubri¬ 
cation regime (Fig. §>)■ 

If one envisages that a second Newtonian plateau ap¬ 
pears at very high shear rates, then a restabilisation 
mechanism emerges naturally. In this situation, the ex¬ 
pected behaviour of the lubrication pressure is shown in 
Fig. §1. Consider increasing the load, starting from a 
small value. At first one follows the right hand branch 
of the pressure curve, until stability is lost at the up¬ 
per saddle-node bifurcation and the system jumps onto 
the left hand branch. Subsequently decreasing the load 
reverses the process, except that the jump back occurs 
at the lower saddle-node bifurcation. The corresponding 
Stribeck curve displays hysteresis and can be constructed 
by noting that the stable fixed points lie on branches of 
the lubrication pressure curve where the flow is Newto¬ 
nian. This is illustrated in Fig. i>. 


B. Normal stress 


Let me now consider some of the other non-linear rhe¬ 
ological effects that might be relevant for non-Newtonian 
lubricants. A long standing question which can be ad¬ 
dressed in the context of the scaling theory is the role of 
normal stress. Again, a simple model capturing this is 


N = 


Vhi 2 

1 + (yr)^ 


( 9 ) 


where ipo ~ t)qt is a normal stress coefficient and /3 is 
another exponent. For concreteness, a typical value of 
this exponent is /3 s=s 1.3-1.5 for polymer melts [ansi- 
The contribution that normal stress makes to the load 
balance can be found by multiplying by the area of 
the conjunction, 7 q ~ Rh 0 . Alternatively we can just 
compare the normal stress with the Reynolds lubrica¬ 
tion pressure. In the low Weissenberg number regime 
(yr < 1) the ratio of these two is 


N 7/> 0 y 2 . h 0 

P VoUro/hl r 0 ' 


( 10 ) 


where the pressure has been taken from Eq. ([2]), and 
U ~ 7/10 and ipo ~ Pot have been substituted in the final 
step. The last factor h 0 /ro ~ (ho/R) 1 ^ 2 is expected to 
be a small number (e. g. < 0.1, say), so the ratio N/p 
is expected to remain small in the yr < 1 regime. In 
particular normal stress should not affect the situation 
in the stable lubrication flow up until the loss of stability 
at the saddle-node bifurcation in Fig. |4j Therefore the 
overall picture shown in Fig. [4] remains unchanged. 

In the high Weissenberg number regime (yr > 1) one 
has iV>Q ~ /zq 1 . As long as /3 > 1 (which is the case for 
polymer melts), this is an increasing function of ho, and 
the gap remains mechanically unstable in this regime. 
Only if /? < 1 will normal stress rescue the conjunction 
from complete mechanical collapse by providing enough 
normal force to support the load in a stable sliding config¬ 
uration. However the above argument shows this should 
not happen until well into the non-linear regime. 


C. Flow transience 

Another factor that can be considered is that the flow 
is transient on a time scale of order the transit time 
t s ~ ro/U. The dimensionless ratio r/t s ~ Ut/vq is 
known as the Deborah number, and its magnitude gives 
an indication of the importance of transient flow effects, 
such as stress overshoot. Substituting U ~ y/io shows 
that the Deborah number is of the order yr x ho/ro■ This 
is the same as the ratio N/p for the normal stress above, 
and the same argument goes through so that the over¬ 
all picture remains unchanged. Like normal stress, it is 
possible that transient flow effects may grow to become 
significant in the non-linear regime since the Deborah 
number scales as h 0 . Hence this may provide another 
mechanism to re-stabilise the gap at a smaller distance. 


D. Extensional viscosity 

In the Appendix, the ratio of the extension to the shear 
components in the lubrication flow is estimated to be of 
the order (h 0 /R) 1 ^ 2 ■ This is the same small number en¬ 
countered above, and here should be compared to the 
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Trouton ratio ( i. e. between extensional and shear vis¬ 
cosity). For example, for polymer melts the Trouton 
ratio is of the order 3 in the low Weissenberg number 
regime (where the gap is stable), and does not increase 
much in the non-linear regime M- On these grounds one 
would not expect extensional viscosity to affect the fric¬ 
tion curve, but the situation may have to be re-evaluated 
if the Trouton ratio is large. 

IV. DISCUSSION 

The main result is that hydrodynamic lubrication by 
non-Newtonian liquids which exhibit strong shear thin¬ 
ning, in the sense that 77 ~ 7 ““ with a > 1 / 2 , is predicted 
to become mechanically unstable as the shear thinning 
regime is entered, with the surfaces closing to near con¬ 
tact. This would be marked by a discontinuous jump 
in the Stribeck (friction) curve as the load is increased 
(Fig. g,). Non-linear effects such as normal stress and 
flow transience are estimated to be subdominant at the 
point of entry into the non-linear regime, so it is unlikely 
they would destroy the basic picture although they may 
provide stabilising mechanisms at high shear rates. Ex¬ 
tensional viscosity is also estimated to be insignificant 
unless the Trouton ratio is large (say, > 10). Depend¬ 
ing on the nature of the restabilising mechanism at small 
gaps, the Stribeck curve may show hysteresis (Fig. §>)■ 

These results were established on the basis of a scaling 
analysis which involves perhaps naive and rather dras¬ 
tic assumptions. For example, the lubrication flow is as¬ 
sumed to remain essentially the same as in the Newtonian 
case, but with a substituted shear-rate dependent viscos¬ 
ity 77 ( 7 ). If this breaks down the analysis would be in¬ 
valid. Furthermore, in Eqs. © and very simple mod¬ 
els were taken for the shear viscosity and normal stress 
as a function of shear rate, and additivity of the nor¬ 
mal stress and lubrication pressure was assumed. There 
are certainly available much more sophisticated (tensor) 
constitutive models, such as the Rolie-Poly equation for 
polymer melts |16j . A major and important avenue for 
further work is to investigate the properties of Reynolds 
lubrication flow using these constitutive models. 

The possible hysteresis in the Stribeck curve makes an 
interesting connection with recent theories of the micro¬ 
scopic origins of shear thickening in non-Brownian sus¬ 
pensions Bang. In these theories the shear thickening 
transition is associated with a ‘breakthrough’ to frictional 
contacts, and indeed the simulations in Fernandez et al. 
m are based on a discontinuous friction curve of exactly 
the form shown in Fig. [4 ]d. Of course, the non-Brownian 
particles that show the phenomenon are more typically 
suspended in a Newtonian solvent, and not a structured 
liquid, and therefore one can question the relevance of the 
present analysis. Nevertheless, experimental systems are 


often prepared with polymer additives (as in Ref. H3, 
and as such they may perhaps exhibit non-linear fric¬ 
tional behaviour of the kind described here if significant 
polymer adsorption occurs. 

I thank M. J. Adams and S. A. Johnson for useful 
discussions. 


Appendix A: Aspects of Reynolds lubrication flow 

Lubrication flow is a superimposition of Couette and 
Poiseuille contributions. For a conjunction between non- 
conformal surfaces, such as that shown in Fig. |T| the 
full solution requires numerics [ 2 ] but basic insights can 
be gained by considering the quasi-one-dimensional case 
M- In that case the flow is 

V ‘ = -\ {h -« ) -Yr i % v(h - V) ' (A1) 

Here x ( cf. r) measures the distance along the gap, y 
measures the distance from the lower surface (moving at 
velocity —17), h is the gap, p is the lubrication pressure, 
and p the viscosity. This velocity field corresponds to a 
net material flux 



This has to be constant even though h varies, and is of 
order UIiq. The implication is that (h^/12p) dp/dx ~ 
Uho, which gives an order of magnitude estimate of the 
pressure gradient. The associated Poiseuille flow velocity 
(along the centerline for instance) is then (h^/p) dp/dx ~ 
U, as utilised in the main discussion. 

The shear rate at the lower surface is 

dv x _ b_ ( A3 ) 

dy y =0 h 2?7 dx 

Since dp/dx ~ pU/h^, the two terms are comparable, 
and the wall stress can be estimated by pU/h 0 . 

The flow field in Eq. © is predominantly in the longi¬ 
tudinal direction. Differentiating with respect to this di¬ 
rection yields an estimate for the extensional component. 
This yields a number of terms, all of which are of the or¬ 
der U/h x dh/dx. Since h = h 0 + x 2 /2R (cf. Eq. Q), 
one has dh/dx = x/R, and setting x ~ Xq ~ (ho HR / 2 
and h ~ ho gives U /xq for the magnitude of the exten¬ 
sional flow component. This agrees with the simple pic¬ 
ture that entrainment involves of the order 100 % exten¬ 
sional strain, in a distance of the order xq, on a time scale 
of the order the transit time xq/XJ. Hence the extensional 
strain rate U/xq ~ ( U/h 0 ) x (ho/xo) ~ 7 x (ho/R) 1 ^ 2 . 
This is the origin of the estimate used in the main text. 
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